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Abstract. 

In this note we present a new propagator for a particle in discrete space under the 

influence of a time-dependent field. With this result we are able to control the shape 
of caustics emerging from a point-like source, as the explicit form of the wavefronts 
can be put in terms of the external field. 
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Similar propagation phenomena in time domain can be found in quantum- 
mechanical systems, electromagnetic waves and sound waves. The analogies between 
their wave equations in controlled and well designed situations have allowed the 
emulation of crystalline structures in settings with periodic symmetry, reaching recently 
a realization of graphene with microwave cavities [T], [2] and photonic crystals [3J. For 
discrete systems without periodic symmetry, the Stark ladder introduced by Wannier 
more than fifty years ago [1] offers itself as an interesting example. In this respect, we 
note that the emulation of electrons under the influence of a constant force has been 
achieved as well: The Wannier-Stark ladder in vibrations of aluminum rods [5] and 
the observation of Bloch oscillations in photonic structures [B] seem to be the simplest 
realizations. Among the most sophisticated, we may single out the propagation of Bose- 
Einstein condensates in periodic optical traps [HI H E] with the possibility of producing 
a Stark ladder by means of a gravitational field. 

In this note we study the more general case of a homogeneous force field modulated 
by an arbitrary time-dependent intensity, with the aim of offering another interesting 
possibility to the already existing configurations and emphasizing the external control of 
the system through such a field. We shall refer to it as a time- dependent Stark ladder, 
pointing out to its generalization through the time dependence of the potential and 
not merely to an equispaced spectrum. Our task is therefore to find the corresponding 
propagator in closed form. In the case of emulations outside of the quantum regime, 
we may simply refer to our result as the Green's function. Armed with the result, 
we shall proceed to characterize the behaviour of caustics emerging from a point-like 
initial condition. We shall also find the explicit relation between the propagation of the 
corresponding wavefronts and the time-dependent modulation of the discrete potential, 
giving the opportunity to discuss the maximal speed of propagation of a signal and how 
it can be controlled within the restrictions of the Lieb- Robinson bound [H]. Two recent 
studies in theoretical [10] and experimental [H] grounds exemplify the relevance of these 
ideas. 

We start with the problem of finding the propagator for a discrete Schrodinger 
equation in the presence of a time-dependent field. One possible approach for 
introducing such an equation is by using the central discretization of derivative 
operators, i.e. 



where a is the lattice spacing, is the mass of the particle and E{t) is the external 
field with dimensions of (energy) /(distance). In solid state physics, one may find a 
similar discrete Schrodinger equation arising from a tight-binding model of a single 
band. Denoting the hopping energy or intersite coupling by A and using the basis of 
atomic functions |0) = Xln*^"!"-) ' "-"^^ 



[0„+i(r) + 0„_i(r) - 20„(r)] + anE(r)0„(r) = ih 



dr 



(1) 



A [0„+i(r) + 0„_i(r)] + E{T)n(t)n{T) = ih 



d(t)n{r) 
dr 



(2) 



Time dependent Stark ladders: Exact propagator and caustic control 



3 



where the field E{t) now has the dimensions of energy. Both (|T| and Q can be simphfied 
by means of a convenient redefinition of units and gauge transformations. In this paper 
we shall work with the tight-binding simplification 

dih (t) 

ipn+i{t) + ^n-i{t) + a(t)n^„(t) = (3) 

which can be obtained from ([T]) through the following definitions: t = — /ir/2/ia^, 
a{t) = 2fxa^E{T)/ri^ and ip^it) = e-^*/Vn(r). 

Some years ago, Yellin [12] found the propagator of ^ for the case a = constant by 
means of the algebraic properties of the hamiltonian. For our problem, the hamiltonian 
operator reads 



H{t)=T + T^ + a{t)N (4) 

with T a discrete translation operator and the position operator with integer 
eigenvalues. The action of these operators on the Hilbert space of atomic (localized) 
functions is given by 

{m^ = ^n+i, (TV)„ = Vn-i, (iV^)„ = n^„ (5) 
and they satisfy the algebra 



[T,Tt]=0, [T,N]=T, [T\N]=-TK (6) 

Now it is evident that for a general function a{t) the hamiltonians at different times do 
not commute: 



[H{tlH{t')] = [a{t)-a{t')][T^-T]. (7) 

In order to obtain the evolution operator corresponding to (|3]), the use of a Dyson series 
seems mandatory. We may circumvent such a cumbersome calculation by writing down 
the Mello-Moshinsky (MM) equations for the discrete representation of the evolution 
operator Un,m (see chapter VII of [13]). Such an operator is related to the discrete 
propagator by the relation Kn,m = 0{t)Un,m- Then, we may obtain K by solving the 
MM equations through elementary techniques for recursion relations. 

It is worth to mention that this procedure was used long ago by Moshinsky 
and Quesne [H] with the purpose to show that linear canonical transformations were 
represented by gaussian kernels. The harmonic oscillator with a time-dependent 
frequency is a good example of this, as it generates a linear canonical evolution and 
its propagator is given by a gaussian in the spatial variables [T5|. Here we employ these 
tools to show that the corresponding discrete version leads quite naturally to Bessel 
functions of field-dependent arguments. 



Time dependent Stark ladders: Exact propagator and caustic control 4 

First, we solve the equations of motion for the operators T,T^ and in the 
Heisenberg picture: 



f = -ia{t)T, = ia{t)T\ N = i{T~T^) 

Using the convenient definitions 



(8) 



f{t)= / rf™(r), F{t)= / rfrexpH/(r)] 
we have the following linear evolution map 



(9) 



T{t) = U^T{0)U = e-*^T(0), 

Tt(t) = f/tTt(0)f/ = e*^Tt(0), 

N{t) = U^N{0)U = N{0) + iFT{0) - iF*T\0). 



(10) 



(12) 



Only (10) and (12) are independent equations, but we include (11) as it shall turn 



useful. We note that A^(0),T(0) are not canonically conjugate opertators, but this is 
not a true obstacle, as they are independent variables and we may employ them in the 
computation of f/n,m- We might equally resort to the canonical pair 



P = ^ (T - Tt) , X = I (T + Tt) iv} 



(13) 



evolving under U, but using (\10h (11 ) and 12) leads to MM equations which are simpler 



to solve. The sought MM equations are the following recurrence relations: 

mm^^n = Um+i,n = e"^^ [t/T(0)]^^„ = e-^^^,„_i 

WmUn = ^Um,n = [U N {Q) + zFt/r(0) - iF*UT\Q)\^^^ 

The first two relations can be easily solved by noting that 



(14) 
(15) 
(16) 



*-^m— l,n— 1 



TT — p-i{n+m)f /2y 
^m,n ^ ' va—ri) 



where V is so far an arbitrary function. Replacing (17) in (16) yields 



(17) 



(m — n)Vrr. 



iFe'f'^Vm-n+1 - iF*e-'f/W^.n-i. 



Finally, we obtain the familiar recursion relation of the Bessel functions [16j by defining 
= arg (F) + 7r/2 + //2, p = |F| and u = m — n. We get 
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P 



K(p)=e''^K+i(p) + e-*<^V;-i(p), 



(19) 



which is solved by K(p) = e~'^'^Jj,(2p), where J is a Bessel function of the first kind. We 
exclude the Bessel function of the second kind due to its irregular behaviour at t = 0, 
violating ?7m,n(0) = 5m,n- Our new propagator reads 



m 



e-''"^Wj„_„(2|F(t)|), 



(20) 



where all the functions can be given explicitly in terms of the external field a, as 
indicated in (|9|. Some properties of (20) can be noted immediately. For example, 
when a = constant, we recover 



sin {at/ 2) 
a/2 ' 



-imf 



(21) 



.\nm 

leading to the usual propagator for the time-independent Stark ladder. Another 
important limit comes from the discrete equation ([T]), where one can let a — > and 
obtain the propagator of a particle in continuous variables. For a detailed derivation, 
see the Appendix. It is also interesting to note that under translations, the propagator 
'picks up' a phase in the form 



Kn+d,7n+d{t) = exp -id / a{T)dT Kn,m{t) 



exp (^-^ ^ {T,Odrd^^ Kn,m{t), 



(22) 



and reminds us of the contribution of the field S to the classical action of a particle in 
continuous space. But the most important feature for the present work is that \F{t)\ is 
not necessarily periodic, opening the possibility of modifying Bloch oscillations by the 
direct control of a{t). 

Our interest now is to characterize the propagation of point-like initial conditions 
in this kind of scenario. Placing the particle initially at the origin gives the probability 
distribution at time t as \Km,o{t)\'^- It is well understood [17J that even in the free 
case, the speed of propagation f in a tight-binding chain has an upper bound (in our 
units fmax = 1) and that the spreading of point-like wave packets occurs at a finite 
velocity. This can be recognized also in our general problem by means of the integral 
representation [IB] of the Bessel function [f] 



|-^m,n(^)| I •An— n (2 1 -f^ I 



1 

2^ 



dk exp [i{n — m)k + i2\F\ sink] .(23) 



I This is not a spectral decomposition of the propagator, in contrast with the free case. 
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From this integral, a fold-type caustic can be extracted. The ray equation obtained by 
the method of stationary phase is 



n — m + 2\F\ cos kn.r 



0. 



(24) 



Solving for kn,m and replacing in the integral of ( 23 ) , gives the curves 

± a/4|FP — {n — my = constant + 2g7r, g G Z. 



[n — m] arccos 



\n — m\ 



2\F\ 



For the main caustic, we simply have zero phase: {n — m) arccos ^ 
■\/4|Fp — [n — mY = 0, with solutions 



|n— m| 



\n-m\ ±2|F(t)| = 0. 



(25) 



± 



(26) 



This brings out the famous light-cones in 1 -|- 1 dimensional space-time and separates 
the propagation (or slow) region from the tunneling (or fast) region. 



ait) 




Figure 1. Left panel: Four possible force fields a{t) producing different types motion; 
the blue curve corresponds to uniform acceleration discussed in example a), the red 
curve corresponds to a mirror placed at p 7^ in example b), the green curve is the field 
for an exponential freeze out of the packet expansion in example c) and the constant 
force in black corresponds to Bloch oscillations. Right panel: The resulting caustics 
from the previous fields. 

Now we turn to the problem of controlling the propagation by means of a. We 
design the shape of our wavefronts in space-time by imposing |i^(t)| and solving for the 
field. This amounts to the inversion of ^ and it is a matter of simple algebra to show 
that 



«W = ^ - -7^. (27) 

p v^- p 

where p{t) = \F{t)\ = ±\n — m\/2 gives directly the position as a function of time. 
The case a = constant can be recovered by setting p as a trigonometric function. 
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Figure 2. Intensity patterns of a point-like distribution placed at the origin at i = 
for different types of external fields. The ordinate represents time and the abscissa 
represents the discrete coordinate. Panel (a): Uniform acceleration. Panel (b): Mirror 
at /5 ^ 0. Panel (c): Freeze out of the propagation. Panel (d): Bloch oscillations 
included as a point of comparison. 



For propagation speeds p > 1, we note that the second term in (27) contains a 



Lorentz factor that becomes imaginary. Turning a into a complex quantity leads to 
non-unitary evolution and exponential decrease, in compliance with the Lieb-Robinson 
bound. In the propagation region, we always have p < 1. Also notable is the presence 



of the acceleration of the wavefront p in (27), which should not be confused with the 



acceleration of a classical particle due to the homogenous field. 

Let us demonstrate the use of our formula with three types of 'engineered' motion 
of the wavefronts: a) uniform acceleration, b) uniform velocity with a perfectly reflecting 
mirror placed at p 7^ and c) exponential 'freeze out' of the wave packet expansion. 

The condition on p for example a) reads p{t) = at'^/2 + vt, where the intial packet 
starts at the origin. This results in a field of the form 



a{t) 



V^l - {at + vy _ a 

aty2 + vt ^1 - {at + vf 



(28) 



The example gives limited motion, since uniform acceleration cannot be sustained forever 
without obtaining an imaginary Lorentz factor. The motion demands an infinite force 
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at times t = 0, {1 — v)/a and a vanishing force at an intermediate time. See the blue 
curves in figjl], where the force and the motion of the wavefront are shown as functions 
of t. The resulting intensity pattern is shown in figj2|a). 

For the motion proposed in example b), we set p{t) = 1 — \vt — 1\ and f < 1. The 
field intensity becomes 



^' l-\vt-l\ ^/T^ ^ ' 

Here it is important to note that in order to achieve a uniform f < 1, we cannot 
simply turn-off the interaction, as this would strictly produce unit speed of propagation. 
Instead, we have a non-trivial solution in (29) even in the absence of the mirror, 
producing a slow but uniform motion. At the event of reflection t = 1/v (mirror at 
p = 1) the field has a delta singularity. For a depiction, see the red curves in figjl] and 
panel (b) of figjS} 

Finally, example c) provides a method to stop the propagation of the pulse by 
applying a field. We take p(t) = 1 — e"'^*, leading to 



= + . (30) 



1 - e-'^* y/l 



The curves for a are shown in green in figjl] and the intensity pattern in figj2]^c). The 
resulting force (30) tends to a non-zero constant as t — ?■ oo. However, such a constant 
field does not produce Bloch oscillations of the packet, since at the origin of time we had 
a very strong (singular) force and the wavefront depends on the history of the applied 
field. One can be convinced of this statement by inspecting the time integrals in (]9|. 

We conclude this note by emphasizing that exact propagators are rather uncommon 
objects [19]. A collection of these kernels has been given in [15] and a number of 
different paths have been devised for their calculation [20] • This has been done mainly 
for continuous problems, including relativistic ones [21]. The discrete case should not 
be an exception. For instance, [T7] contains a Feynman path integral version of the 
free discrete kernel. What we have presented here is a method that is common to both 
continuous and discrete realms and that explains the resulting solvability for a wide 
class of systems: The spatial representations of canonical transformations. 

In a less technical order of ideas, let us mention that the examples presented here 
are of an illustrative character. However, we should not disregard completely their 
applicability to experiments designed ex professo, which may range from the simplicity of 
torsional waves in a piece of metal to the ambitious control of Bose- Einstein condensates 
in a dilute regime. 
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Appendix: Continuous limit 

Here we compute the continuous limit of the propagator corresponding to the 
Schrodinger equation ([T]). Let us start by writing our kernel in the appropriate 
units: According to the definitions, we have r = —2jjLa^t/h, tpn = e*'*^/'"" 0„ and 
a{t) = -2ma^E{T)/h^. This leads to 



Kfirn (^) 



where 



\F{t) 



exp 



-i ( m/(r) + 



J^_„ {2\F{t)\) 



(Ai; 



/(r) = ^ dsE{s) 



Fir) 



n 



2fia- jQ 

The continuous limit corresponds to 



ds' exp 



dsE{s) 



(A.2) 
(A.3) 



a — )■ 0, an X, am — > x', dx = a, (A. 4) 

while other quantities such as fi, t and E{t) remain fixed. For our computations we 
need the following ascending expansions in a: 



Fir) 



hr 



2/ia 



/ ds / dsE(s) + - — - / ds' 
Jo Jo 4:fxh Jo 



/ ds' 


[ dsE{s) 


Jo 


Jo 



-\ 2 



(A.5) 



\F(t) 



hr 1 

+ 



2/xa2 Ajjhr 



-\ 2 



ds' / dsE{s) 
Jo 



Ajjh Jq 



ds' 



1 2 



dsE{s) 



arg[F(r)] ^ --^ Tds' T dsE{s). 
"''T Jo Jo 



(A.6) 



(A.7) 



By substituting (A.6) in (A.l), we observe that the limit a — )■ demands the use of an 
asymptotic form of Jn{z)- Such an approximation comes from the Meissel expansion 
^16j and has the form 
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Jniz) 



V2 



TTIZ 



exp 



i[z + 



2z' 



(A.8) 



which corresponds to 1 << n << z. With the expansions (A. 5 - A. 7) and the asymptotic 



form (A.8), we finally get the limit 



K{x, x\ t, 0) = dx 



X exp 



X exp 



27rihT 

ix' 
'T 



exp 



2hT 

dsE{s) - 



exp 



2fih Jq 

PT PS' 

I ds I dsE(s) 
Jo 



i{x' — X 




hr 






ds' 


If" 




Jo 



ds' / dsE{s) 
Jo 



-\ 2 



dsE{s) 



(A.9) 



In passing, we note that the resulting asymmetry in x, x' is an effect due to time 
irreversibility of the external field; the propagator necessarily depends on the two 



variables x+x' and x — x'. One can verify that (A.9 ) is the correct limit of the propagator 
by checking that it satisfies the equations of motion of the continuous problem. In fact, 
the MM equations yield 





OX ox' 



ihr d 



K{x,x';t,0) 
K{x,x';t,0) 



- / dsE{s) 



- f ds' f dsE{s 
At Jo Jo 



K{x,x';t,0) 

K{x,x';t,0) 



(A.IO) 
(A.ll) 



/i dx 

which are equivalent to the following Heisenberg equations of motion 

P{t) =Po- dsE{s) 
Jo 

T 1 r , r' 

x{t) = xq -\ — po / ds' dsE{s). 

f^Jo Jo 

Furthermore, we can recover the well-known result for a particle in a constant 
homogeneous field by letting E{s) = Eq. We have 



(A.12) 
(A.13) 



lElr^ 



24/i/i 2^hT 



n 2 



ds' / dsE{s) 



2fih Jq 



ds' 



dsE{s) 



i{x + x')EoT i{x' -x) r , , r , ^ ix' r 



(A. 14) 
(A.15) 



and with this result, the phases in (A.9) reduce to the ones reported in [15], formula 
(6.2.18), page 175. 
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